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Abstract 

Quark-hadron duality in the inclusive semileptonic decay B XJv in the Shifman-Voloshin limit 
A <C Sm = mil — rric <C rrib, rUc is studied within a nonrelativistic potential model. The integrated 
semileptonic decay rate is calculated in two ways: first, by constructing the Operator Product 
Expansion, and second by a direct summation of the exclusive channels. Sum rules (Bjorken, 
Voloshin, etc.) for the potential model are derived, providing a possibility to compare the two 
representations for T{B Xclv). An explicit difference between them referred to as duality- 
violation effect is found. The origin of this effect is related to higher charm resonances which are 
kinematically forbidden in the decay process but are nevertheless picked up by the OPE. 

Within the considered order the OPE and the sum over exclusive channels match each other, 

up to the contributions of higher resonances, by virtue of the sum rules. In particular this is true 
for the terms of order Srri^ /rri^ and ASm/m^ which are present in each of the decay channels and 
cancel in the sum of these channels due to the Bjorken and Voloshin sum rules, respectively. The 
size of the duality violation effects is estimated to be of the order 0{K^^^ /m'j.dm^) with 6 > 
depending on the details of the potential. Constraints for a better accuracy are discussed. 
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Quark-hadron duality in the inclusive semileptonic decay B Xdv in the Shifman-Voloshin limit 
A <C 5m = mi, — mc <C mi,,mc is studied within a nonrelativistic potential model. The integrated 
semileptonic decay rate is calculated in two ways: first, by constructing the Operator Product 
Expansion, and second by a direct summation of the exclusive channels. Sum rules (Bjorken, 
Voloshin, etc.) for the potential model are derived, providing a possibility to compare the two 
representations for T{B —f Xclv). An explicit difference between them referred to as duality- 
violation effect is found. The origin of this effect is related to higher charm resonances which are 
kinematically forbidden in the decay process but are nevertheless picked up by the OPE. 

Within the considered order the OPE and the sum over exclusive channels match each 

other, up to the contributions of higher resonances, by virtue of the sum rules. In particular this is 
true for the terms of order 6rr? Irrii and t^5m/m^ which are present in each of the decay channels 
and cancel in the sum of these channels due to the Bjorken and Voloshin sum rules, respectively. 
The size of the duality violation effects is estimated to be of the order 0{K^^^ /m^5m^) with 6 > 
depending on the details of the potential. Constraints for a better accuracy are discussed. 

I. INTRODUCTION 

The interest in inclusive decays of heavy mesons is two-fold: experimental study of such decays can provide im- 
portant information on the weak mixing angles of heavy mesons, and a theoretical treatment of such processes which 
includes also nonperturbative effects is possible. The theoretical framework based on combining the Operator Product 
Expansion (OPE) and Heavy Quark (HQ) expansion provides decay rates and differential distributions as series in 
inverse powers of the heavy quark mass with the coefficients proportional to the matrix elements of the operators of a 
proper dimension A remarkable property of this expansion is that in the leading-order this is just the free-quark 

decay, and the first correction appears only at order l/mg. 

On the other hand it is understood that the quark-hadron duality technically implemented through the OPE is an 
approximate framework ||]. For example, the calculation based on OPE does not take into account all details of the 
hadron spectrum which lead to the dependence of the set of open decay channels on the momentum transfer. The 
OPE ignores this fact and this inevitably yields some errors in the OPE results 

The theoretical description based on the OPE represents the decay rate as a contour integral in the complex g^-plane 
(for details see the next section). The OPE can be justified only in regions of the complex g^-plane away from the 
physical region, whereas in the case of the calculation of the decay rate (both differential and integrated) the contour 
always involves a segment which is close to the physical region |^]. This can lead to duality- violating effects, i.e. the 
difference between the exact and the OPE based results. 

However, it is not easy to estimate the errors arising in the OPE, since the exact hadron spectrum in QCD is 
complicated and not exactly known. So, testing directly the accuracy of the quark-hadron duality is only possible in 
few exceptional cases. Examples discussed in the literature are QCD in the Shifman-Voloshin (SV) limit and the 
't Hooft model §. 

In the 't Hooft model (2-dimensional QCD with Nc — > oo) the spectrum is reduced to an infinite number of single 
bound states and known precisely so that the direct summation of exclusive channels is possible. First numerical 
analysis of the sum over exclusive channels reported the presence of the duality- violating 1/toq correction for the 
total width . Later the summation was performed analytically for the case of a massless light quark . The result 
of the OPE calculation agreed with the exact result in this case through 1/toq order. 

Duality in QCD in the SV limit |^] has been studied in I pd] , ^^ . This limit requires Aqcd ^ 5m — mi, — iric ^ toq . 
A peculiar feature of the SV limit is that a summation over exclusive channels becomes possible due to kinematical 
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reasons: the process occurs near the zero recoil and thus only few decay channels contribute in the leading 1/toq 
order. The expansion of the relevant transition form factors in this kinematical region is known and the sum over 
exclusive channels can be evaluated. The absence of Aqcd/"1q corrections to the free-quark result in the semileptonic 
(SL) decay rate has been demonstrated in |llT|] . However, to check the absence of Aqcd/'tiq and Sm/mQ corrections 
within the SV kinematics is not enough to ensure duality in the 1/toq order in the general case, beyond the SV 

limit. Namely, one should also check that potentially large terms of order O ^Aqcd'^™"/"^q^^) which are present 
in individual decay rates cancel in the sum over exclusive channels. The analysis of the Kqcy^Sitl / ttiq terms in the 
exclusive sum was performed in |]l2| for QCD in the V — A case. It was found that the duality within this accuracy 
requires a new sum rule. The full comparison to higher orders has not yet been performed. 

We study the quark-hadron duality in the SV limit within a nonrelativistic potential model. The model has several 
features which make it especially suitable for this purpose: the model is self-consistent in the SV limit; the spectrum 
of bound states is relatively simple and can be calculated; the exact representations of the transition form factors in 
terms of the hadron wave functions are known. These features provide a possibility to calculate the exclusive sum. 
We adopt a technical simplification of a Lorentz scalar current instead of the V — A current, like it is done in ref. 

The main purpose of our analysis is to check whether or not the OPE result calculated to some order is equal to 
the sum over exclusive channels expanded to the same order. Both series are double expansions in powers of A/we 
and A/ 6m. They are asymptotic series and the question of their convergence is left for a later publication |l3| . 

Our main results are as follows: 

• We construct the expansion of the T-product of the two currents in a series of local operators (the OPE) in 
the potential model for a general form of the quark potential. Technically this is done by the expansion of the 
Lippmann-Schwinger equation. We consider the expansion to all orders in A/ dm but neglect terms ~ A"/m" 
with n > 3. This OPE series provides the expansion of the differential and integrated semileptonic decay rates 
in powers of A/mc and A/ 5m. 

Let us point out that the OPE series in the potential model has an important distinction from the Wilsonian 
scheme in the field theory: Namely, in QCD (perturbative) contributions of small distances below the scale is 
referred to the Wilson coefficients while contributions of large distances above this scale is referred to the matrix 
elements of the local operators. As a result both the c-number Wilson coefficients and the matrix elements of the 
local operators acquire the /z-dependence. In the potential model we also expand the average of the T-product of 
the two current operators over the B meson in a series of local operators, but the resulting c-number coefficents 
as well as the average values of the local operators (see Eq.(||)) do not have a scale dependence. 

• The OPE and the sum over exclusive channels are related to each other by sum rules, similar to the Bjorken 
[p^ , the Voloshin and the whole tower of higher moments We derive these sum rules. They involve an 
infinite sum of terms corresponding to all hadronic excitations, with each term having a well-defined heavy mass 
expansion. The question of the heavy mass expansion of the sum (in other words, of the uniform convergence of 
the series) has not been tackled in this paper. If the contribution of higher excitations vanishes rapidly enough 
the uniform convergence is expected. 

• The OPE provides a heavy mass expansion for the inclusive semileptonic decay rate. To compare it with 
the result of summation over exclusive semileptonic decay channels we make use of the sum rules. An explicit 
difference between the two expressions is found, both for the integrated and the differential rates. This difference 
corresponds to the contribution of the resonances kinematically forbidden in the decay process which are picked 
up by the OPE. This 'unphysical' contribution is related to the poles in the complex q°-plane outside the physical 
region which however contribute to the OPE result. The size of this duality-violation cannot be estimated in 
all generality since it depends on the potential and on the convergence properties of the sums over resonances. 

• For the integrated decay rate the OPE prediction and the sum over exclusive channels match, up to the duality- 
violating contributions of higher resonances, within the order: Terms of order Sm^/mq, (5mA/m^, which 
are present in any individual decay rate cancel in the sum over all channels thanks to the Bjorken p4| and 
Voloshin |^5| sum rules, respectively. For terms of order A^/mg, A"^ /mgdm, etc the agreement {again up to 
contributions of higher resonances) is provided by the higher moment sum rules. The duality-violation induced 
by the kinematical truncation of these higher resonances in general has the order 0(A^+''/m^(5m'') where b 
depends on the details of the potential V{r) both at large and small r. 

For the smeared differential distributions near maximal the violation of the local duality is found at the order 
ASm/ml. 
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We make an explicit proof of the present results for the special case of the harmonic oscillator potential in Ref . . 
This is important since some demonstrations given below are rather formal. 

In the next Section we present some details of the kinematics and discuss the analytical properties of the decay 
amplitude. In Section 3 the 1/toq expansion of the quark propagator is performed and the OPE series for the SL 
decay rate in nonrelativistic quantum mechanics is constructed. In Section 4 we consider the HQ expansion of the 
exclusive form factors in the potential model, and derive the inclusive sum rules (Bjorkcn, Voloshin, etc.) which are 
crucial for comparing the exclusive sum and the OPE result. In Section 5 we provide an analytic expression for the 
duality- violation contribution and identify its origin. We estimate the accuracy of the OPE both for the integrated rate 
and the smeared distribution near zero recoil. A special emphasis is laid on discussing the role of different inclusive 
sum rules in establishing the relationship between the OPE and the sum over the exclusive channels. Conclusion 
summarizes our results. 

II. KINEMATICS AND THE ANALYTICAL PROPERTIES OF THE DECAY AMPLITUDE 

We consider the inclusive SL decay B XJi^. The rate of this process reads 

r{B ^ = ^eiq' > mmwiq) (1) 



2Ms J 2tt 

where L is the leptonic tensor, and the hadronic tensor W is defined as follows 

d'px0{p°x)S{Px - M],){B\J\X{p)){X{p)\J+\B)diipB - Px - q)- (2) 

X 

Here the relativistic normalization of states is implied 

{p\p')^2p"i27r)H{p-ff). (3) 

For the sake of clarity we assume the technical simplification that the leptons are coupled to hadrons through the 
scalar curren^ In this case the leptonic tensor is a scalar function of only one variable, g^, and the hadronic tensor 
W depends on the two invariant variables v = Pb ■ q/Mg and q^. In the rest frame of the i?-nieson these are g*^ and 
q^. At (7° > and fixed q^ the sum in (||) runs over the hadronic states with masses Mx < Mb — \/q^- The decay 
rate can be written as follows 

r{B ^ XM) = ^J dq^dq°\mq° > mL{q^)Wiq°, q^), (4) 
2Mb 

with g2 = (g0)2 _ -2_ 

Equivalently, we can use q° and q^. Let us consider the W{q°,q^) as an analytical function of at fixed q^. One 
can write the following relation 

ij^w{q°, f) = -Ldisc,or(g°, f) (5) 



where 



T{q\f) = I dxexp-*«-(i3|r(J(x), J+(0))|S) = ^r^ffi^^^, (6) 

' 2Mb J \ I V w: \ ))\ I 2Mb ^ Mx - Ex{-q) - q° 



Ex{—q) = y/q"^ + M'^ is the energy of the state with the mass Mx and the total 3-momcntum —q. The sum over 
X in Eq. (0) for T runs over all hadron states with the appropriate quantum numbers. The selection of the states 
kinematically allowed in the decay process is made by the proper choice of the integration contour in the complex 



^ Recall that for the case of the V — A current and massless leptons, the leptonic tensor has the form L^,^ ~ dfivQ^ ~ 9^9" 
and for the scalar current L ^ . We consider throughout the paper the leptonic tensor of the generalized form L — {q^)^ . 
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plane. Namely, the decay rate r{B XJv) can be represented as the contour integral in the complex plane over 
the contour C{q^) which depends on the value of (Fig 1) as follows 

r{B ^ XJl.) = ^. J dq'\q] J dq^e{q''>\<^\)L{{qy-q^)T{q'>,e)- (7) 

C(g-2) 



c 



/ 










e — e — Qy' 




h3 — e — e — e — e — e — e — ¥b 

^- «S ' 





FIG. 1. Singularities of the amplitude T{q^,q^) in the complex q'^ plane. Circles are hadronic (eg) poles which are only 
singularities in the confined potential model, and the cross stands for the free b c quark process. The vertical line Re (g") = |g| 
separates the kinematically allowed region of the real axis from the kinematically forbidden region. The contour C{q^) embraces 
all states in the allowed region. Poles at the left of the boundary correspond to kinematically forbidden bound states. 



It is important that the left crossing of the contour C{(f') with the real axis in the complex q^^ plane occurs at the 
point 90 = I?!; otherwise the contour can be freely deformed in the region where the function To(goj9^) is analytic. 
We symbolically mark this constraint with a 0-function in the integrand. The integration over such contour selects 
at any given q"^ only physical states which can be produced in the decay B XJv, i.e. states with the invariant 
masses such that A/|. + q^ < Mb — \q\- Notice that whereas the left crossing of the contour with the real axis is 
tightly fixed at the point Re (q°) = Im (g^) = 0, the right crossing of the contour with the real axis can be safely 
moved to the right. In the general relativistic case there are cuts which correspond to other physical processes. In 
the SV limit these cuts are separated from the physical decay cut by windows of the width O^ttiq). In the potential 
quark model such cuts are absent. 

The amplitude T{q'~', q^) is given by the T-product of the two local current operators, which is the classical case for 
performing the OPE. Namely, one has 

Tiq°,e} = ^TT^ / rfa;cxp(-zgx)(i?|T(J(x), J+(0))|i?) =^C„(g°,g2)^5|0^j5^^ (g) 
2Mbi J ^ 

where 0„ are local operators and Cn{q^,q^) are the c-number coefficients. Introducing the expansion (||) into (|^) 
gives the integrated rate as a sum over various local operators. 

We shall obtain the integrated SL decay rate within our model by the two means: first, we construct the OPE 
series for T{q'^, g^), and second, we calculate directly the sum over exclusive channels. 



III. THE MODEL 

We consider this decay in the SV limit, 

A <^ 5m — nib — "m-c <C "m-c, "rrib- (9) 

Notice that in a non-relativistic model, A refers to a fixed energy scale proportional to the light quark mass m^, 
to the average quark momentum in the hadron rest frame ((i3|A;^|i3))^/^, to the parameter (3 defined in Eq. ( [2^ ) 
and hbd or es to be defined in Eq. (|ll|). These parameters may be strongly hierarchized, for example a genuine 
non-relativistic situation implies {{B\k'^\B)Y/'^ <C m^, but all these quantities remain constant as mc, mf, oo, they 
remain proportional to some fixed hadronic scale which we call A by analogy with QCD. This is to be distinguished 
from 5m which is taken as an independent parameter. Thus we consider the double limit 5m/mc — > 0, and A/6m — s- 0. 
Notice finally that \q\ is of order 5m. 

To avoid confusion, it is important to stress that the standard OPE expansion assumes 5m /mq constant, even if 
small. So the order of a term 0(A"+™/m"((5m)™) in this paper corresponds to the order 0(A"+™/mg"'"'") of the 
standard OPE expansion. 
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We treat the leptonic part relativistically, but for the calculation of the hadronic tensor we use the nonrelativistic 
potential model. The nonrelativistic treatment of the hadronic tensor is consistent within the SV kinematics and can 
be used as a tool for studying some of the aspects of quark-hadron duality. We shall make use the fact that in the 
nonrelativistic potential model we know the structure of the hadron spectrum and have an exact representation for 
the hadronic matrix elements of the quark currents. 

It is convenient to use the nonrelativistic normalization of states (which is used hereafter) 

{p\p')^{2nf5ip-j/), (10) 

and consider the process in the rest frame of the decaying i?-meson. The B meson is the ground eigenstate of the 
Hamiltonian Hm, 

Hm\B) = Mb\B) = (mfc + jud + eB)\B). (11) 
In the _B-rest frame this Hamiltonian has the form 

Hbd = mi, + md + h h Vbd{r) = nib + ma + hd- (12) 

Zirib 2.rad 

For the B Xc transition we need the cd bound states with the total 3-momentum —q, which we denote Dn{—q). 
These are eigenstates of the Hamiltonian 

ffcd(g)=m, + mrf + ^^-t^ + -^ + I4rf(r). (13) 
Zrric ^nid 

such that 

HcdmDn{-q)) = EDM\Dn{-q))- (14) 

In this equation Ed^ {q) is the nonrelativistic energy of the bound state Dn with the 3-momentum —q 

Ed^ (q) = Md„ + —, ; r , Md„ = nic + md + eD^- (15) 

2(mc + md) 

The expression (H) for the decay amplitude T now takes the form 

T(,o,,-)^El^4.-)l %,^_^^(^-^_^o ^ (16) 

where Fn{q) is the B —>■ £)„ transition form factor, 

Fn{q) = {B\J\Dn{-q)), Mb ^ nib + ma + es, 



and the sum runs over all cd resonances. The expression (16) can be also written as 

T(gO,g) = {B\J- ^ , J^\B), (17) 

Mb - g° - H,diq) 

where {Hcd{q) — E)~^ = Gcd{q,E) is the full off-energy-shell Green function (propagator) of the cd system. The B 
decay amplitude is thus given by an average of the Green function Gcdiq^E) at the point E = Mb — q^ over the B 
meson. 

Let us specify the transition current operator Jb^c- For the sake of argument we neglect the quark spin effect and 
consider spinless nonrelativistic quarks and choose the quark current in the form 

Jb^c = / dkdk'b{k')c+{k), (18) 
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where c(6) is the annihilation operator of the c{b) quark 0. 

For the quark current (^8|) the B — > _D„ transition form factor in the rest frame of the i?-meson reads 

Fn{q) = f dkg iJB{kg)^D^ (kg + q) , (19) 

where kg is the momentum of the hght spectator. 

Similarly, for the current jl^ ) the expression (|l^) takes the form 

^('^"''^^ = <^It7 ,^ 1-^)- (20) 

Mb - - H,d{q) 



IV. THE OPE OF THE DECAY RATE 

The main idea in constructing the OPE series for T, Eq.(^, is to single out Hm from Hcd{q) in the denominator 
in Eq. (|2^) and to use the eigenvalue equation (pi]). First, let us introduce the operator 6H{q) which measures the 
difference of the denominator of Eq. (E(J) from the inverse Green function of the free-quark transition 



Explicitly, one finds 



Mb-q"- H,M = {5m-^- 9") - 5H{q). (21) 

Irric 



SH{q) = Hcd{q)-mc~md- es- (22) 

Next, isolating hi,d in SH{q) we obtain 

6H{q-) = ih,d -eB) + U—- —) {P + V,) + ^+0 (^) , (23) 
2 \mc nib J rric \ m% J 

where the scale (3 is provided by the hadronic matrix elements 

P"" ^{B\k''\B)^{B\Vi\B). (24) 

As already mentioned, (3 is of the order of A. The quantity Vi here is a part of the expansion of the potential Vqq in 
powers of l/mq 

Eqs. (pO|) and (|2^) allow us to construct the expansion of T{q'^, q"^) in inverse powers of 6m — g° as follows 

Making use of Eq. (ESh we obtain 



■^Notice that the standard scalar current reads Jb-,c~J ■^^b{k')c^ (k) and in the nonrelativistic limit takes the form 
~ / dkc+{k)dk'b{k') {l - {l - . Neglecting the factor (l - (l - as done in Eq. (|l|) leads to technical 

simplifications both in the OPE and in the exclusive sum. It can be easily realized that a particular choice of the current 
however does not touch any arguments related to duality. 
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T{q^,(f) = {B\B) 



1 



2m c ^ 



(B|(/iM-eB)|B) 



1 



5m — if- q'^ 



(B\ 



\B) 



3 oo 
i j'^l n— 



( 5m — TT 

\ 2771 a ^ 

IB). 



3+n 



+ 0(AVm^). 



(27) 



The remainder has the order 0(A^/m^) if we keep ~ Sm^ and go fixed. 

Finally, using Eq. (|l]) and the relations {B\k,\B) = and {B\k,kj\B) = \5ij{B\k'^\B) we find the following OPE 
series 



I 



Sm-4^-qO 2ml 



[B\k^ + Vi\B)- 



1 



27nc 



■E- 

n=2 



{B\On-2\B) 



T + 0(A7m^), 



(28) 



where 6„ = Yl^j=i kj{hbd — eB)"kj. Hereafter the J2 symbol is omitted. We denote {B\k^ + Vi\B) = Pq, (3o ~ A. 

The series (|2^ ) is a double expansion of T{q^, q^) in K/nic and A/ {5ra — q^), limited to second order in A/rric 

and expanded to all orders in A/ {6m — (f). The poles are at 



qi[r)^^^- 



2mr 



(29) 



The first term in (j^^) gives the free quark decay amplitude. A remarkable feature of this series is that the A/ 5m and 
A/toc corrections to the free-quark decay are absent thanks to Eq. ( pi] ) and the relation {B\ki\B) = 0. The expansion 
( Pq ) substitutes the whole set of hadron poles by a complicated quark singularity at the point q^ = q^{q'^). 

Let us treat the series ( |2^ ) formally and calculate the integrated rate which is obtained as a double expansion in 
A/mc and A/6m. 

Let us rewrite the expression (Bh) as follows 



.0 



T{q",q 



1 - i^A!±il^.™- ' 



2ml 



dSm 



{B\0„.2\B) 1 
3m? n! 



1 



f (5m - - (70 



(30) 



This representation is very convenient for the calculation of the decay rate Eq. (|^): the integration over q° is now 
easily performed since 



dq"L{iqy-q^)9{q^>\q\) 



C{q^) 



(dm — 4— 't 

\ 2mc ^ 



L{ql{f), f)e (\(f\ < -me + ^ml + 2mc(5m) , (31) 



where the ^?-function 6 [\q\ < —mc + ^ m^ + 2mc(5m^ reflects the fact that the left crossing of the contour with the 
real axis in the complex q'^ plane should always happen at the point Re((/*') = |(z|. The integrated rate is then given 



by the expressior 



^ A remark is in order here. When computing the decay rate in Eq.(^^ we have interchanged the derivation with respect to 
Sm and the integration over dq'^. We can also directly integrate the expression Eq. ( psj ) over lig". In this case we should take 
into account that the contour C(g^) in the complex g'^-plane al way s lies on the r.h.s. of the line Re(g°) — \q]. If we erroneously 
do not take this condition into account then multiple poles in (Bq) do not contribute at all since the complex integral vanishes 
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where 



_/ {B\k' + V^\B) 
2m? 



^-^ 3m^ n\ \ dSm J 

n—2 *- ^ 



(Sm.mc), (32) 



For the free quark decay one finds 



df\q\^L{[q",if)]^~f 



r(fe — > civ) — Ii[5m,mc). 



(33) 



(34) 



Let us consider the leptonic tensor of the general form L{q^) = (q^)^ . For semileptonic decays to massless spin 1/2 
leptons N — 1. The case = corresponds to scalar leptons. Since the leptonic tensor is proportional to q^, it is 
now convenient to introduce a new integration variable q^ as follows 



Then the integrated rate takes the form 

T(^pE(^B ^ XJi^) = r(6 ^ du) 

{B\k^ + Vi\B) 



(35) 



(36) 



2m2 



dSr 



-Ii{Sm,mc) + ^ 



{B\d^^2\B) If d Y 



n=2 



3m? n\ V dS' 



h{5m, ruc), 



where we have taken into account that Ii{Sm,mc) gives the exact free-quark decay rate. A simple algebraic exercise 
gives to the 1 /m^ accuracy 



Ii{Sm,mc) = (<5m)2^+3 
/3((5m,mc) = ((5m)2^+5 



^f.o (1 - ^ 



3 Sm 15 Sm^ 



H/2 



+ o 



2mc 8 m? y 8 ^'^ m, 
Sm 
mr 



where 



A. 



N 



f dxx-"(l - = i?(iV + 1, m + 1), = ^]^^V2' 



i3(p, g) being the Euler function. Finally, we come to the relation 



T{B ^ XM _ ^ I {B\k^\B) ^ {B\Vi\B) 



T{h civ) 



2m? 



2m? 



n— 3 



2A^ + 5 mzdrri 



(37) 



(38) 



(39) 



with C„ = k\(n-k)\ ■ Notice that the coefficient of the term {B\k \B) does not depend on N, i.e. it does not depend 
on the form of the leptonic tensor. 

Summing up, the OPE predicts the following features of the inclusive SL decay rate: 

• The LO term reproduces the rate of the free-quark decay process h ^ c. 



when the multiple poles are inside the contour as well as when they are outside (see also discussion in ]2|). To proceed correctly, 
one should replace a multiple pole, say a double pole, by an equivalent set of two neighbouring poles. Then the crossing of the 
boarder gives a non-vanishing result when one of the two poles is inside the contour and the other one is outside. Integrating 
before taking the derivative with respect to 5m as in Eq. (^) corresponds to treating in a specific way the crossing of the 
boarder by the multiple poles. We show elsewhere that both treatments lead to the same result. 
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• The l/jTic and l/5m. corrections arc absent. This is due to the fact that the average over the S-state of the 

operator of the relevant dimension vanishes. 

• Lowest-order corrections to the free-quark process emerge in the l/m^ order. A main part of these corrections 
is due to the average values of the dimension-2 operators {B\k'^\B) and {B\k'^ + V\\B). Also the operators 
On = kj{hbd — eB)"'kj contribute in the 1/m^ order. Their contribution is however suppressed with the additional 

powers of Sm. 

In the next section we shall analyze the accuracy of the OPE predictions. 



V. HEAVY QUARK EXPANSION AND THE HADRONIC SUM RULES 

Before proceeding with the direct summation of the exclusive channels one by one we derive hadronic sum rules 
which are important for the comparison of the exact result with the OPE analysis. 



A. Heavy quark expansion of the form factors in the potential model 
The wave function of the Qq bound state has the form 

rnQTUq (kg kq]] 



rriQ+mq \ nig tuq 



^pikq^kq) = S{p~kq- kq)'tp 

= S{p- kq - kq)^P { kq ) (40) 

V mq +mq J 

where p is the momentum of the bound state. 

The B Dn transition form factor is the average over the meson states of the operator ^bc{(l) given by the 
following kernel 

{h\^bcmkc)=5{h-h-q). (41) 
So the transition form factor is defined by the following expression 

{B{pB)\^bcmDn{Pn)) = S{PB - Pn - q)Fn {{vb - Vnf) , (42) 

with Vb = Pb/ {rrib + mg) and w„ = p„/ (mc + m,) and 

rriq _ \ , It' rUg 



Fn (vb -Vn)= I dkqtpB kg ; PB i^D^ kq ; Pn ■ (43) 

J \ rrib + mg J \ rric + mg J 

A simple change of variables kq kg + mq/{mh + mq)pB makes it obvious that the decay form factor depends on the 
square of the relative 3- velocities of the initial and final mesons, and not on the relative 3-momentum squared as the 
elastic form factor. Nevertheless in the B rest frame we write 



Fn{q) = J dkgi)B{kq)'4>Dn {k, 



mc + rfiq 



(44) 



The wave function i/'Qg is an eigenstate of the Hamiltonian 



^Q9-=^fz^ + :^)+^Q«W, (45) 



2 \mq m, 



where k is the conjugate variable to r. The transition form factors Fn has some general properties independent of 
the details of the potential Vqq. Such properties of the transition form factors are derived by performing the HQ 
expansion of the Hamiltonian. To this end we apply the usual Quantum Mechanical perturbation theory. 

For our purposes it is convenient to consider hbd as the full hamiltonian, hcd as a nonperturbed Hamiltonian, and 
U = hbd — hcd as the perturbation. The perturbation has the form 
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where we assume the foUowing expansion of the Qq potential 

Vqq = 1^0 + T^Vi + -^V2 + ... (47) 
2mQ 2mg 

The perturbation has the order Sm/ml such that we can construct the HQ expansion of the wave functions and 
the binding energies. Let us remind the standard formulas: Let {V'-D„} be the full system of eigenstates of the 
hcd, and the {eD„} the corresponding eigenvalues. Then, the mass of the rt-th excitation in the cd system reads 
Md„ — rric + nid + e_D„. Let {ipB„} be the full system of eigenstates of the hbd, and the {es^} the corresponding 
eigenvalues. 

The standard formulas give 

= ^D,. + V + 0{5m^pymt) (48) 

^ ^D^ — eD„, 

and 



= + C/„„ + V + ... (49) 

where Q 



f/m„ = -TT^i^AJ^' + l^il'AD^) +0 (50) 
2mj \ m-^ ) 

The excitation energies satisfy the relation 

e_D„ - eu„ ~ (n - ■m)Xnm, >^mn ^ /3. (51) 

In terms of the wave functions, the transition form factor (Q) takes a simple form 

Fniq) ^ ii^Bol^DA-q}) (52) 

The expansion of the wave function ipBo = "0-8 reads 

V'Bo = V-Do + V o + 0{dm^Pyml), (53) 



such that 



^ 2 \mb rricj (Do - f-D™ 



/on + V /mn C) + 0((5m /3 /mj 54) 



where jmn[(i) = {''pD^\'>pD^{<t)) ■ By virtue of (^ one obtains 



/nn('Z) = 1 - rlr,^ + 0(g"'/m^), 

fnmio)^ rlm—^ + 0{q^/m'^), m^n, (55) 



*eso is just es defined in Eq. jl We also usually write B instead of Bq. 
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with rmn being numbers of order unity plus higher order 1 /rric corrections. We shall use the notation r„ = r„o ■ Notice 
that the radii r„ describe the form factors of the transitions between different levels in the cd system (Dq — s- Z?„) and 
so know nothing about 5m. 

We now rewrite Eq. (^J) as follows 

At (f ^ = we thus come to the relation 

Fo{0) = 1 + 0{Sm^(3^/mj). (57) 

One can see that the 0{5m[3 /m^) term in i^„(0) is absent. This is a non-relativistic analog of the Luke theorem [Q. 
For the squares of the form factors we obtain the following important relations^ 

FS{q)^l~pl^ + 0{—^), pl = rl+0l'" 



\ J \ ITT-] 

F^iq)^ pI^ + 0[ 1^ , Pl^rl + Oi^). 58 

As we shall see later the radii p„ (as well as r^) are not independent and satisfy certain sum rules. The relations (|5^ 
are the main result of this section. They are necessary for the calculation of the decay rates. 

B. Inclusive hadronic sum rules 

To obtain a non-relativistic equivalent of the whole tower of sum rules [p^ , i.e. the Bjorken sum rule, the Voloshin 
and the higher moments, we consider the following set of quantities (i = 0, 1, . . .) 

SM = {B\i5Hiq)y\B), (59) 



where 5H{q) is defined in Eq. (21). Notice that Si{q) appear in the expansion for T{q^ , (f), Eq. (26). We shall derive 
two different representations for Si{q) and obtain sum rules equating these representations. 

The first representation is obtained by inserting the full system of the eigenstates \Dn{—q)) of the Hamiltonian 
Hcd{q) in Eq. (|59| ). The |_D„(— g)) are also eigenstates of the operator 5H{q) that is made obvious using 5Hcd{q) in 
the form Eq. (^21): 

5H{q)\D,,{-q))=5,mDn{-q)), J„(g) ^ ezj„ - £b + ^^^J^ ma) ~ 
As a result of inserting the full system we find 

oo 

5.(q) = 5]|K(q)r(J„(g))\ (61) 
ri=l 

Eq. ( ^ ) gives the following expansion for 5n{(t)'- 

Sni^f) = A„ + 1 - ^) iDo\k^ + y.l^o) - , , + O ( ^) , (62) 

2 \mc rrib J 2mc[mc + nid) \ m'^ J 

where 

A„ EE eD_^ - (63) 



^ At any n states with angular momenta L = 0, ...,n exist. The form factors fn{<l^) and F„{q'^) are thus understood as 
properly normalized sums X]l=o-/'"^® ^^'^ '^l^o^nLiq), respectively. 
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Notice that within the leading-order accuracy we can replace {Do\k'^ + Vi\Dn) with {B\k'^ + Vi\B). 
Another representation for Si{q) is obtained by using SH{q) in the form (p3|): 

SM-{B\{SH{q)y\B)^{B\(h,-es + ^^(- -^)+tl+o(^)] \B). (64) 

\ 2 \mc mbj nric \ m-^ J I 

This formula gives Si(q) in terms of the matrix elements of various operators over the B- meson. 

The representations ( |6l|) and ( pi]) for Si{q) provide the l.h.s. and the r.h.s. of the sum rules, respectively. Let us 
notice that terms denoted by 0(/P(5m/m^) in Eqs. ( p2[ ) and ( p^ ) do not depend on q. All q'-dependent terms are 
shown explicitly. 

Using Eq. ( pT| ) and the relations {B\ki\B) — and (B\kikj\B) — ^Sij{B\k'^\B) we come to the set of sum rules. In 
fact each of these sum rules is equivalent to an infinite number of relations at different powers of and l/rric- 
i=0: 

OO 

So^Y.\Pn{q)\' = l. (65) 

n 

Obviously the r.h.s. does not depend on q. At q^ — this is an identity. Using the definition ( p8| ) and comparing the 
term linear in q^ we find the NR Bjorken sum rule 

OO 

Pl = J2pn- (66) 

n=l 

i=l: The r.h.s. of this sum rule reads 

S. = U^~^)iB\k^ + V.\B)+0(^). (67) 



2 yrric 

where we have used Eqs. ( pT| ) and (|l^ . The r.h.s. of this sum rule is also independent of q. From the definition (|5^ 
and using (|6^) as well as the SR (|65|), we rewrite the Eq. ( |67| ) as follows 

E . ^ . + O (^) . (68) 

Notice that the terms O{0^5m/m'j.) cancel between r.h.s amd l.h.s. Comparing the linear in q"^ term yields the NR 
Voloshin sum rule 

HpI^^-^YT^- (69) 

n=l mc 

Let us notice that the r.h.s. of the Eq. (|6^) does not contain higher-order l/rric corrections. 

Combining the Bjorken and the Voloshin sum rules provides a simple constraint on the parameter which is in 
fact the slope of the Isgur-Wise function. Namely, 



■2 



n— 1 n— 1 n— 1 



where A„ are defined in Eq. ( |63[ ) 

i=2: The r.h.s. of this sum rule reads 



S...fiEip^ + 0(^}. ,71) 



3m 

Using (H) and (|^) yields for the l.h.s. 

q^{B\P\B) f ^f(3^STn\\ q^ml ^ f P^Sm 

n=0 



Efi^(g1A^-^^W^ 1 + + , 7 ■ (72) 

— ' 3m^ \ \ J J Am-^(mc + md) ^ ' 
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The linear g^-term yields 

f:Pl^l-liB\k'\B)(l + o(^-^)) (73) 



i>3: 

For i=3 we find for the r.h.s. 



.3 = imMhti^MM^ (i+o(^))+o ) . (74) 



3 V \ m 

Using (|5|)-(0) yields for the l.h.s 



3m? V V m? 

n— \ \ L 

-4 /^2fl2A^\ -6 /«5j^2 



The linear g^-term yields 

'|fcj(fefcd 

3 V V "^r 7/3 V V m^ 



E P^A3 = (^|fc.-(^M-e.)fc,|i?) + o f = i(B|Or|B) f 1 + O f ^ ) ) . (76) 



n=l 

t2 



Similarly at higher i > 3 one obtains at the pSm/m^ accuracy 

f^plK = \{B\k,{hM-eBr~%\B) = i(i?|4-2|i?). (77) 

n—l 

These sum rules are used in the next section for comparison of the exact decay rate with the OPE result and for 
analyzing the duality-violation effects. 



VI. SUMMATION OVER THE EXCLUSIVE CHANNELS 



We now proceed to the summation of the exclusive channels. As the first step, let us show that there is an explicit 
difference between the exclusive sum and the OPE series. 



A. The origin of duality violation 

Proceeding with the sum over the exclusive channels we write: 

2TTt ^ J Jc(q^) Mb-q^ - En{q) 
= ^E / 'i^'^i'^') I dq^\'^l jf-^^T 

= 7r-2^ dqL{q) dq \q\- 1-- + ... 

= ^ fdq^Liq^iy: f [i 'Jl^ + .1 , (78) 

where in the r.h.s. \q\ = a/ (9")^ — q^- Notice that the sum is truncated at the proper n{q^) which is the maximal 
number of hadron resonances kinematically allowed at a given value of g^, i.e. resonances satisfying the relation 
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Mn < Mb — The contour C{q^) is responsible for this selection, since only the resonances enclosed by the 

contour contribute into the sum. All states which are beyond this contour do not contribute. 
Finally, the series (78) can be written in the form 

r r "^^'^ 

r{B ^ XM) = J dq^L{q')e{q') J dq^dq^miq^ f - <f - f) E 



Siq" 



2mr 



6m). 



On the other hand, the sum rules (67 )-([<' 4|) allow us to rewrite the decay rate (|3^) in the form 
rOP^^(B^X,M= J dq'Liq')0{q') J dqVm{<ir - - 



dSm 



1 9^ 

2 oom'^ 



6iq" 



2mr 



6m). 



(79) 



(80) 



It is easy to see that the exact result and the result of the OPE are different due to contributions of highly-excited 
states: at any q^ the OPE picks up also the contribution of the resonances forbidden kinematically at this q^. Thus 
the accuracy of duality is determined by the accuracy of violating the sum rules connected with the truncation of the 
exclusive sum, and is therefore connected with the convergence of these sums. 



B. Sum of the exclusive channels and the accuracy of the OPE 

We now calculate the individual decay rates keeping terms of order ^ (^)" the decay rates but neglecting 
higher orders ^ (^)". The necessary expressions with the relevant accuracy are given below. 
\q\ in the free-quark decay b —^ civ at q^ has the form 



I ^ y/^ — 2 2 (i _|_ ^ ^^^^ ^ 9^ \ 

\ 2mc 8 TO^ 8™^ / 



(81) 



The general expression for \q\ in the B Dnlv transition at q reads 



JM„ 3 6Ml 6M^-q^ 



|<zl„ = ^6Ml-q^ 1 - „. " . + " , + ^. " (82) 

2[mc + md) 8[mc + m^ 8(mc + md)^' 



where (5M„ = Mb - M„ 6m - A„ - 6mPo/{2mc) from Eqs. and (|e 
The necessary accuracy for the transition into the ground state is 



/ , , / P?, \ , f 6m 3 6m^ 6m^~q^\ 

Recall that /3g = {k^ + Vi). 

For |(j|„^o less accuracy is enough since the contribution of the -D„,n 7^ into the SL decay rate is suppressed by 
the additional factor g^/m^: 



\q]„^0 = \/i6m-Ar,)^-q^. (84) 

With these formulas for \q] the decay rates of the exclusive channels take the following form. 
Free quark decay b civ: 

1 dT{b civ) 5m ^ 3 5m'^ q^ \^ 

2 c A^2 _ Jl 



L{q'^) dq^ \ TOc ' 2 m^ 2m1 J 



rz — ^ ^ / i 5m 15 Am^ b 6m^ — q \ 



14 



The B — > DqIi/ channel 



using the definition (|5f 

The S ^ DJv {n ^ 0) channeh 



TO 

Now everything is ready for the calculation of the integrated SL decay rate. We again consider L = {q^)^ 



1. The integrated rate and the global duality 

It is convenient to represent the results for the partial decay rates in terms of their ratios to the free quark decay 
rate. The latter has the form 



r(6 ^ civ) = (,5to)2^+3 
Making use of the relation (JSq) we find 



I 3 6m 15 5m^\ 5 ^ (5to^ n f 



T{B ^ Dolv) ^ _ 3pg Sm^ 3 rud 5m _ /3g 
T{b^clv) 2N + 5 ml 2 1 + md/mc ml ^ '2ml 

r{b cliy) ml 2iV + 5 V <5m . 



(89) 



k=2 



Some remarks are in order: 



1. The main part of the OPE (i.e. the free quark decay) is reproduced by the T{B DqIv), within the leading 
and the subleading I/toc orders. The excited states contribute only within the {Sm)'^/ml and ASm/ml orders 
in the SV limit. 

2. Nevertheless, each of the individual exclusive channels contains potentially large terms of the order Sm^ /ml 
and ASm/ml which are absent in the OPE series. 

Now summing over all exclusive channels we find 

r{B XJiy) ^ _ 5n^ f 3 _ "y^^" 2^ 3 q ( 1 "^rf _ x ■ -z . 
T{h-^clv) ml I 2N + 5 ml [21 + md/mc 



2ml 2ml 



1 i~i)'c',^^, {3j::=t pia^^) 

mlSm ^ 2A^ + 5 8m^-^ ^ ' 

k—Z 

The sum over the charm resonance levels is truncated at n„^ax j which is the number of the resonance levels opened 
at q^ = 0. For the confining potential and in the SV limit rimax is found from the relation A„^^^ ~ 5m. 

Using the sum rules (|6^)-(77|) to rewrite the OPE result ( ^9| ) as the sum over hadronic resonances, the difference 
between the OPE and the exclusive sum (the duality-violating contribution) explicitly reads 

r(b ^ civ) = 2N + 5 + 
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where 



J2 Pli^nf^ E [rl + 0{K5m/ml)\{^^f. (93) 

As expected, this duahty-violating contribution is connected with the charm resonance states forbidden kinematicaUy 
in the decay process. This kinematical truncation of the higher resonances induces a violation of duahty equal to 

i^^^^jt*^) for every fc < 2iV + 5. 

To estimate the error induced by the truncation and thus the size of the duality-violation effects, we need to know 
the behavior of the excitation energies and the transition radii at large n. 

1 . For quite a general form of the confining potential we can write the following relations for A„ for large n (recall 
that in the SV limit A„^^^ ~ 5m) 

A„ > ACn", n > Umax, (94) 

with C and a some positive numbers. In particular, this estimate is valid for the confining potentials with a power 
behavior at large r. 

This estimate for A„ is only depending on the behaviour of the potential at large distances (the infrared region). 

2. The transition radii satisfy sum rules similar to sum rules for in Section V, namely^ 



oo 



md 



2(1 + md/mc)' 

oo 

Er2(A„)('=+2) = -(i?o|A:,(/icd-eDo)'%|i?o), A; = 0, 1, 2, . . . . (95) 

n=l 

Hence, the behavior of the radii at large n are connected with the finiteness of the r.h.s. of the sum rules. We can 
guarantee this for the Bjorken and Voloshin sum rules, where finite values stand in the r.h.s. (the ground state radius 
rg is finite for the confining potential). In general, the finiteness of the matrix elements of the operators kjihcd — 
(such as e.g. the kinetic energy of quarks in the ground state) depend on the properties of the potential at small r 
[the ultraviolat behavior) and probably also at large r [the infrared behavior) Q We have assumed throughout the 
paper that the average kinetic energy of the light spectator quark in the ground state is finite, i.e. (-Do|fc^|^o) — A^. 
This already restricts some properties of the potential at small r and provides convergency of one more sum rule in 
Eq. (|95|). If, in addition to this, we assume that the average values of the operators kj{hcd — e_Do)'^fcj for k — 1, ...,K 
over the ground state are finite, then combining with the behavior of the energies at large n we come to the following 
estimate 



1 



This allows us to obtain the duality- violation originating from the truncation of the various sum rules: 



e > 0. (96) 



Bjorken : ^^(") . ^ y < ^ M ^ ^ ^ f _^ ^ " 

Voloshm: ^ V < A -— ^-^ 7- ■ (97) 



^ Notice that these relations are exact and do not have any l/rric corrections. 

'^We don't have a classical Wilsonian scheme where the ultraviolet region is referred to the Wilson coefficients and the infrared 
region is referred to the matrix elements of the operators, so we can have these regions mixed 
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Similar estimates can be done for higher moment sum rules. One can see that the truncation error in any of the sum 
rules leads to the duality-violation of the same order O (A^+^/to^Jto^) . An interesting feature about these estimates 
is that the dependence on a has disappeared from the final result. Hence, the estimates are independent of the details 
of the potential at large r, provided the potential guarantees the confinement, i.e. a is positive. 

These are however rather crude estimates which do not take into account further possible suppressions (due e.g. to 
the orthogonality of the wave functions of the ground n = and the excited states n > 0). In such a case the real 
accuracy is better, and might depend on the details of the potential also at large r. In general, we can state that the 
truncation (duality-violaiton) error occurs at the order 



m 



8m) ' 



where the exponent 6 > depends on the properties of the potential (in general, both at short and long distances). 
A more detailed analysis of which potentials satisfy the above requirements is beyond the scope of this paper and is 
left for another publication |p^ . 

If we would like to have the truncation error of a higher order in 1/mj,, e.g. in ©(A'^/mj?), this is not so straight. 
Namely, in this case we need 

0[t-t-\- (99) 



8ra ySm'^mc 

As we have noticed, the series for (J'-''-' in the main part does not depend on rric, so the only possibility to have the 
relation (^9|) fulfilled in the SV limit, is to have = starting from some number n. (Exactly this situation takes 
place for the HO potential where all Dq Dn transition radii for n > 1 are equal zero [|7j). In this case for large 
enough values of 5m, the term proportional to in ( p3|) disappears and the second term provides the truncation 
error of the order 0{Sm^A^/m'^). 

As we are going to show elsewhere, the accuracy of duality of order A^/m^ can be achieved if we keep a fixed ratio 
5m/ mc when mc — > oo. One can proceed exactly along the same lines, but technically a bit different treatment is 
necessary: namely, at several places throughout the paper terms of the order 5m^ /m^ have been omitted, and they 
should be kept if the limit 6m/ mc — const is considered. This analysis will be presented in p3[ |. 

Finally, it is interesting to notice that all resonance levels opened at = are contributing on equal footing to 
the sum rules and therefore to the decay rate. So, a considerable delay in opening channels with large n compared to 
the channels with small n with the increasing recoil does not matter at all. This is a very irnportant feature which 
basically determines a high accuracy of the OPE calculation of the integrated decay rate (cf. [l6|Jl2|). 

2. The smeared q'^ -distribution and the local duality 

The situation however differs considerably if we consider the differential decay widths. We find it more physical 
to use here the four vector variable. The region near q'^^ix (zero recoil) is special: as we move to higher g^, the 
excited channels close one after another leaving ultimately only the Dq ground state opened. 

Let us consider a partially integrated decay rate in the g^-region above the threshold of the Dn=i channel. In this 
case the relation between the OPE and the exact result (which is reduced in this case to the exclusive B DqIv 
decay) reads: 



dq^ - J "^^^ 

((5m-Ai)2 (5m-Ai)2 



1 + 



Ai(5m 



(100) 



In this formula we have negl ected a difference between the upper boundaries of the quark and hadron channels of 



the order A^/m^. Eq. (100) means that local duality near maximal q^ is violated at order 0{Ai5m/m^). As we 
have seen, the dangerous terms of this order are cancelled in the integrated rate against similar contributions of other 
channels due to the Voloshin sum rule. However, the 0(A(5m/m^) violation of the local duality might have negative 
consequences for the application of the method to the analysis of the experimental results. For example this happens 
if one observes only a small part of the phase space near maximal M. 



17 



VII. CONCLUSION 



We have studied quark-hadron duality in decays of heavy mesons in the SV limit using the nonrelativistic potential 
model for the description of mesons as qq bound states. Our main results are as follows: 

1. The OPE is constructed and the following A/rric and A/ 6m double series is found for the integrated decay rate: 

where Cfe are calculable constants and fco depends on the leptonic tensor. 

2. The HQ expansion of the transition form factors in the nonrelativistic potential model is performed. A non- 
relativistic analog of the Luke theorem for the exclusive transition form factor between the ground states is 
obtained. 

It is shown that the sum of the squares of the B Dnliy transition form factors are expressed through the 
expectation values of the operators emerging in the OPE series. These nonrelativistic analogs of the Bjorken, 
Voloshin, and higher order sum rules provide a bridge between the sum over exclusive channels and the OPE 

series. 

3. The integrated decay rate is calculated by direct summation of exclusive channels. For the comparison of this 
directly calculated r{B — > XJi/) and the corresponding T'~'^^{B — > X^ly) the sum rules are necessary. A 

difference (duality-violation) between the two expressions is observed. As shown explicitly by the use of the 
sum rules, this difference is connected with the higher cd resonances which are forbidden kinematically in the 
decay process but are implicitly taken into account in the OPE approach. Therefore the accuracy of the OPE 
is directly related to the error induced by the kinematical truncation in the sum rules (Bjorken, Voloshin, etc). 
The actual error depends on the convergence of the series, i.e. on the nature of the potential. We have discussed 
the constraints on the latter convergence which lead to the duality violation of order 0{A'^~^^ /rrt^Smf') with b 
depending on the behavior of the potential both at the short and long distances. 

Up to the mentioned duality-violation, the agreement between the OPE and the exclusive sum is achieved within 
different l/rric orders due to different reasons: 

The leading order and the subleading Sm/rric and A/rric orders the free quark integrated decay rate T{b civ) is 
equal to the rate of the transition into the ground state D^. This is due to the specific behavior of the transition 
form factor between the ground states near the zero recoil (Luke theorem). Also part of the A^/toq correction 

in the OPE result proportional to the {B\l? + V\\B) matches to the contribution of the ground state Dq in the 
exclusive sum. 

For higher order terms the agreement between the OPE and the exclusive sum is a collective effect due to subtle 
cancellations in the sum over exclusive channels: 

Namely, each of the individual decay rates V{B Dj,lv) contain potentially large terms of the order 5rn? /m^ 
and A5m/m1. These terms cancel in the exclusive sum due to the Bjorken and Voloshin sum rules, respectively. 
The higher order sum rules allow us to represent the contribution of exclusive channels in terms of the average 
values of the operators Oj over the B-meson state. 

4. If the differential semileptonic decay widths are considered near maximum q^, the violation of the local duality 

occurs at order 0{A5m/m\). 

Clearly, in QCD the situation is more complicated because of the multiparticle states, pion emission, hybrid 
and multiquark exotic D mesons, radiative corrections. Nevertheless the duality violation due to the kinematical 
truncation of the series should be quite similar to the case of nonrelativistic quantum mechanics. Also similar is the 
role of the inclusive sum rules in obtaining the duality relations. 
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